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ABSTRACT 

Thin non-thermal X-ray filaments are often seen in young supernova remnants. We 
used data from the 1 Ms Chandra observation of Cassiopeia A to study spectral prop- 
erties of some of the filaments in this remnant. For all the cases that we examined, the 
X-ray spectrum across the filaments hardens, at about 10% level, going outward, while 
observed filament widths depend only weakly on the photon energy. Using a model that 
includes radiative cooling, advection and diffusion of accelerated particles behind the 
shock, we estimated the magnetic field, turbulence level, and shock obliquity. 

Subject headings: shock: diffusion — advection, supernova remnant 



Introduction 



Young s upernova remnan ts (SNRs) have long been thought to be the main source of galactic 
cosmic rays (jShklovskii Ill953l ). Evidence for the existence of high-energy electron s in SNRs first 



came with the detection of non-thermal emi ssion in the radio and later in X-rays (e.g.. lKovama et al 
19951 : bamba et allbood : Isiane et allboOlh . With the us e of the Chandra X-Ray Observatory, de- 
tailed images of SNRs have revealed very thin structures (jLong et al.ll2003l ) near the forward shock. 
The spectral and spatial properties of such structures, often referred to as filaments, are consis- 
tent with synchrotron emission from highly relativistic electrons. High-energy protons and nuclei 
are also believed to be produced within SNRs, although no direct evidence has been conclusively 
found. However, a recent analysis of the Fermi-LAT spectrum of the SNR W51C suggests that the 
main component of e mission in the GeV band from this object is produced through interactions of 
high-energy hadrons (lAbdo et al.ll2009l ). 



In SNRs, charged particles may gain 
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off m agnetic turbulences downstream and magnetic irregularities upstream (jBlandford fc Eichler 
19871 ). Such diffusive shock acceleration (DSA) naturally leads to a power-law distribution of 
particles, which is seen for cosmic rays. However, little is known about the precise nature of 
the magnetic turbulences. The scattering of particles by circularly polariz ed MHD waves with 



a fre quency equal to the gyrofrequency of the particles results in diffusion (jBlandford Eichler 



19871 ). Quantifying the magnetic diffusion, therefore, offers a way to probe the turbulent wave 
spectrum in SNRs. 



Many issues regarding DSA are still unresolved (for a summary, see iRevnolds 1 120081 . and 
references therein). For instance, some of the observed electron spectral distributions have indices 
that are either too high or too low to be accounted for by the mechanism of DSA alone and may 
require that other effects be considered, such as n onlinear proce sses (for example, the deceleration 
of the incoming fluid with respect to t he shock (jEichler Ill979l ) or the excitation of MHD waves 
by the accele rated particles themselves ( Bell 1978a! ) ). as well as the effects of second-order Fermi 
acceleration (jOstrowski fc Schlickeiser 1 11993), 



The injection problem is another unresolved issue. One of the requirements of DSA is that 
particles have enough initial energy to be able to pass through the shock without being significantly 
deflected and, therefore, that their gyroradii be greater than the shock thickness. Since the electron 
gyroradii in the thermal plasma of SNRs are typically smaller than the shock thickness, this poses 
a problem for ex plaining the initial accel eration mechanism for these particles (for a discussion on 



this problem, see lMalkov &; Drury 



2001 



The understanding of D SA requires knowledge of magnetic field. There is evidence for magnetic 
field amplification in SNRs (jBerezhko et alJbood : Ivink k Laming] I2OO3I : feerezhko & Volk Il2004h 
with respect to pre-existing ambient field. A mechanism for explaining such amplification has 
been treated for the case where the magnetic field is parallel to the shock normal (i.e., with an 



obliquity angle of zero, see lBell fc Lucek 11200 ll ). although much is still unknown about the possible 
dependence of this amplification with obliquity angle, shock speed, or composition of the upstream 
medium. 

Another feature of shock acceleration that has been debated is the so-called shock precursor, 
a region supposedly formed by scattering of high-energy particles upstream of the shock. Since 
DSA requires that particles cross the shock back and forth it is possible that particles with very 
high energy scatter in a region ahead of the shock. It is expected that the scale of the precursor 
will depend on properties such as the particle density and the diffusion co efficient. The exist ence 
of such a precursor is still debated and its properties are largely unknown ([Ellison et al.lll994l ). 



To gain insights into some of these unresolved issues we studied X-ray synch rotron emissio n 
from Cassiopeia A (Cas A), a young SNR with an age of approximately 300 yr (jHughes Ill980l ). 
Much has alre ady been learned ab out this remnant, from the identification of the forward and 



reverse shocks ( Gotthelf et al 



2001 



to the measurements of the mag netic field (|Vink fc La 



2003) as well as the observed proper motion of the forward shock (IDeLaney &: Rudnick 



mmg 



2003 
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DeLanev et al.l 12004 ; iPatnaude &: Fesen 1 12009 ) . The expansion rate of the SNR seems to be lo wer 
than expected, due probably to more efficient particle accelera tion (IPatnaude fc Fesen 1 120091 ) or 
to a more complicated density profile of the ambient medium ([Hwang Laming 1 12009 ) . As will 
be seen, our analysis does not allow us to draw conclusions regarding these possibilities, partly 
because it only deals with the leptonic component of cosmic rays in Cas A, and partly because we 
do not consider the dynamics of the expansion. 



In this work, we used data from the 1 Ms exposure on Cas A with Chandra (jHwang et al.l l2004) 
to carry out a detailed spectral analysis of non-thermal filaments in the outermost region of the 
remnant. Specifically, we are interested in the energy dependence of the width of the filaments and 
spectral variation across them. In the context of advection and synchrotron radiative cooling, the 
widths are expected to decrease with increasing energy and the spectrum softens going downstream. 
However, diffusion may significantly modify the behaviors. 



The Observation 



We extracted spectra from nine non-thermal filaments of Cas A from the archival Chandra 1 
Ms observation, taken with the backside-illuminated S3 CCD chip of the Advanced CCD Imag- 
ing Spectrometer (ACIS). The observation consists of nine segments, the first was taken in 2004 
February and the rest between 2004 April and May (see iHwang et al.l 120041 ). For our analysis 
we used level 2 archival data products and reduced the data with the standard Chandra software 
package CIAO v3.4 and Chandra calibration database (CALDB) version 3.5.2. Data were accumu- 
lated in GRADED mode to avoid telemetry loss, therefore the effects of charge transfer inefficiency 
in the spectra cannot be corrected. This loss of charge affects the measured pulse-height distri- 
bution and the energy resolution, although the effect is small for the backside-illuminated CCDs 
(|Townslev et alJboonh . 

The image of Cas A, shown in Figure [H was obtained by combining events in the energy range 
from 0.3 — 10 keV and then correcting the count map by effective exposures. Since the effective area 
is energy dependent, weighted exposure maps were calculated at different energ ies and combined . 
We focused on regions that had previously been identified as being non-thermal (|Stage et al.l l2006) 
and been thought to be associated with the forward shock (jGotthelf et al.ll200ll ). However, due 
to low statistics (even with a 1 Ms exposure), faint filaments mainly located in the western and 
eastern sides of the remnant were not included in the analysis. The image also shows the nine 
filaments chosen for this work along with off-source regions for background estimates. 



2.1. Dividing the filaments 



Each of the non-thermal filaments was divided into an 'inner' and an 'outer' region, with the 
'inner' region being closer to the interior of the remnant. The division between the inner and outer 
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regions for each filament was set at the peak of its linear intensity profile. This division made it 
possible for us to quantify the difference in the spectral properties of the radiation emitted by the 
electrons at different locations. One may naively attribute the difference to the fact that electrons 
in the inner region have had more time to evolve after interacting with the forward shock than the 
electrons in the outer region. However, as we will show, there seems to be a fair amount of mixing, 
implied by the inferred diffusion coefficients. 

As an example, Figure [2] shows the linear profile of Filament 5 (for the 0.3 — 10 keV band). 
This profile was obtained from a 1".5 x 6". 9 region running perpendicular across the filament with 
a bin size of 0".5. The top panel shows the division between the inner and outer emitting zones. 



2.2. Energy dependence of filament widths 

The widths of each filament were estimated in three energy bands: 0.3 — 2 keV, 3 — 6 keV and 
6 — 10 keV. In order to quantify the width of a linear profile, we fitted the profile around the peak 
with a Gaussian function, as shown in Figure El The results of the fits are summarized in Table [TJ 
No strong dependence of the widths on energy is apparent. 

We should note that the overall linear profiles of the filaments are highly non-Gaussian (see, 
e.g., Figure [2]). Nevertheless, we think that the derived Gaussian widths reflect fairly accurately 
the widths of the filaments. 



2.3. Spectral evolution of filaments 

To carry out the spectral analysis, we reprojected the event 2 files to a common tangent point 
and used the CIAO tool acisspec to extract events between 0.3 keV and 10 keV from each region 
shown in Figure [H calculate weighted Auxiliary Response Files, and combine the spectra from the 
individual segments of the observation. Consistent calibration was used separately to produce the 
exposure-weighted responses by applying time-dependent gain corrections appropriate for — 120°C 
GRADED mode data on the back-illuminated S3 chip. After the individual files were combined, we 
binned each spectrum such that each bin c ontained at lea st 100 counts and proceeded to individually 



model them with XSPEC version 11.3.2 (lArnaud 111999 ) 



All spectra show, with varying degrees of prominence, the presence of emission lines (see 
Figure indicating the existence of thermal photons in the regions. For most cases, we added 
two Gaussian components to model the lines at around 1.85 and 2.38 keV, which we attribute to Si 
K XIII and S K XV, respectively. Other weaker lines also appear to be present in some filaments. 
Filaments 5 and 6 show two additional lines at 1.3 keV and 1.0 keV, most likely associated with 
Mg XI and Fe XXI, respectively. Filament 9 also shows the line at 1.0 keV. 

The inner and outer extraction regions have a typical extension of about 7" each, which cor- 
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responds to a phy sical size of roughly 0.1 pc (assuming a distance of 3.4 kpc to the remnant; (see 



Reed et al. Ill995l )). The regions for Filaments 2 and 5 are smaller (4" and 3", respectively) since 
there seems to be a considerable amount of thermal emission in these areas. To assess possible 
"contamination" from thermal emission, we also experimented with thinner extraction regions for 
each filament as well as on-source background regions. In the first case, we failed to remove the 
lines seen, while in the second one it becomes difficult to determine the appropriate locations of 
background regions necessary to avoid subtraction of non-thermal photons. The resulting lack of 
statistics after the subtraction generally does not allow to carry out a satisfactory analysis of the 
non-thermal X-rays. It is possible that the thermal and non-thermal emissions are cospatial, but 
we do not rule out that the detection of thermal photons might be due (at least partly) to scattered 
X-rays. 

The photon spectra of all filaments were satisfactorily fitted with an absorbed power-law, with 
indices ranging from 2.2 to 3. Figure |4] shows the spectral fits, as well as residuals, for both the 
inner and outer regions of each filament. The results are summarized in Table [2j 

Although in most cases the error intervals of the photon indices for the inner and outer regions 
overlap, we note that the spectrum of the inner region is in general softer than that of the outer 
one. The difference in photon indices between the inner and outer regions is on the order of 10%. 

The hydrogen absorption column values obtained from the fits are typically 0.7—0.9 x 10 22 cm~ 2 
in all regions except for a larger value of 1.3 x 10 22 cm -2 found for Filament 8, at t he western edge 



of th e remnant, where it is believed that it is interacting with a molecular cloud (jKeohane et al 



19961 1 



3. Theoretical modeling 



We developed a model to explain the observational results. The model takes into account 
synchrotron radiative losses and diffusion of particles in the forward shock region. We assumed 
that the injected particles follow a power-law spectral distribution with index T and the particle 
spectrum subsequently evolves. 



We approximated Cas A as a sphere with radius R = 10 19 cm (jReed et al. II1995I ). The non 



thermal emission is assumed to come from a thin shell near the edge of the sphere and integrated 
along the line of sight. The evolution of the non-thermal electron distribution is given by the 
diffusion-loss equation. We used the solution derived by Syrovatskii (1959), but also included the 
advection process. For Cas A, the advection speed of the plasma downstream of the shock is 



V, 



adv 



1300 km s~l, equal to the shock speed divided by a shock compression ratio of 4, 
(^sh ~ 5200 km s~ 1; Vink et al. 1998). This value agrees with X-r ay Doppler shift measur ements, 
which imply a velocity relative to the shock of about 1400 km s~l (jWillingale et al.l 12002). 



We note that we neglected energy loss due to adiabatic expansion, because it is expected to 
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have little effect on the distribution of particles in thin filaments. Similarly, energy losses due to 
Bremsstrahlung radiation and inverse Compton are not considered, since the synchrotron loss is 
expected to dominate. To derive analytical solutions, we approximated the synchrotron radiative 
power for an electron as 

P»{l) = (cttcB 2 ^ /6tt)5(v - u e ), (1) 

where v c = 3qB~/ 2 / 'Airmc = 3z^7 2 /2, and vl is the Larmor frequency; here, 7 is the Lorentz factor 
of the particle, m is its mass, q is its charge, and B is the magnitude of the magnetic field. We will 
discuss the effects of this approximation in Section 5. 

We assumed that all of the emission originates behind the shock, where the magnetic field is 
believed to be amplified with respect to the ambient field. The diffusion coefficient was taken as 

where k is a proportionality constant to be determined. The case when k = 1 is referred to as Bohm 
diffusion. O ther types of diffus i on are also being studied, including Kolmogorov and Kraichnan 



turbulences (jKolmogorov I Il94ll ; iKraichnan I Il965l ). but will be discussed in detail elsewhere (D 



Lomiashvili et al. 2010, in preparation). 

There are four main parameters in the model: magnetic field, spectral index of electrons (r), 
diffusion length (tajf) and advection length (^ a( j v )- The diffusion and advection lengths are defined 
as 

/ 3 \ 1/2 
/ Kmc \ 

w = w 

respectively, where tp = aT/Q^mc. For convenience, we combined these quantities to define two 
new parameters, Ajj£ = l^/R and ( = i a ^ v (^-keV)/l gj ^ v , which can be determined from the data. 



4. Results 

We implemented the model in XSPEC as a table model and applied it to the spectra of the 
filaments. For each spectral fit, the line features and the hydrogen column density were fixed to 
values found in the corresponding power-law fits. The quality of the model fits is the same as that 
of the fits by this phenomenological (power-law) model. 
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4.1. Magnetic field and diffusion coefficient 

From the best-fit A^y? and £, we derived the magnetic field and diffusion coefficient for each 
filament (Equations (3) and (4)): 

B „ 56,0 (^Y 2 " (±Y 2/3 ( % r V /3 (5) 

P 10.02/ V 5.0/ V 1.3 x 10 8 cms- 1 J \ 10 19 cm ' 



« « 0.05 — sov . (6) 

V5.0 / V 1 - 3 x 10 cms" 1 / v ; 

The results are summarized in Table El The indices of the injected electron spectrum found vary 
from 2.6 to 4 and the magnetic field ranges from ~ 30 (iG to 70 fiG, while k stays around 0.02 
(but can reach up to 0.1 — 0.15). We should stress, however, that the error bars for the last two 
quantities are considerable. The magnetic fields for Filaments 1 and 7 are most uncertain, due to 
the difficulty in constraining (. 



4.2. Maximum energy of electrons 

The accelerated electrons will lose their energy due to synchrotron radiation. The evolution of 
the particle's Lorentz factor, 7, is given by 

- d i) =-^V (7) 

A maximum energy will be reached by the particle when this loss becomes equal to the accel- 
eration rate. If we assume that the mean magnetic field is perpend icular to the s hock normal, then 



for a compression ratio of 4 we can write the acceleration rate as (jJokipii 1119870 

i^i = (8) 
7 dt J acc 32kD b ' V ; 

where Dg = r g c/3 is the Bohm diffusion coefficient and r g = (mc 2 /qB)j is the particle gyroradius. 

Our assumption about the direction of the magnetic field is justified by our estimates of the 
diffusion coefficient, which constrain the obliquity angle to be nearly 90° (see below). 

The maximum energy for an electron then is given by 

E — - < 660W > (03s) ~" 2 (30^) " 2 ( 5 .2xitfw0 ■ (9) 
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4.3. Shock obliquity and turbulence level 



We considered diffusion in the radial dire ction, with a corresponding diffusion coefficient given 



by (jJokipii 



1987 



Blandford fc Eichler Ill987h : 



D = D\\ cos 2 6 + D\ sin 2 



(10) 



where 9 is the angle between the mean magnetic field and the normal direction of the shock. Here, 
we assumed that the kinetic theory relations, Dn = ijDb for the diffusion coefficient along the 
mean direction of the field and D± = t]Db/(1 + rj 1 ) for t he component of th e diffusion coefficient 
perpendicular to the mean direction of the field, hold (e.g.. lForman et al.lll974l ). In these equations, 
r) = A m fp/r 9 is the particle's gyrofactor, which is the ratio of the scattering mean free path, A m fp, 
to the particle gyroradius, r g (jHavakawa Hl969i ; lMelrose 1 11980). Since isotropic Bohm-type diffusion 
is assumed here, we can rewrite Equation (2) in the form D = kDb- From Equation (10) we have 

/ 9 „ sin 2 9 \ „ . 

rj ( cos 2 6 + rr—^ ) ■ (11) 



When diffusion is taken as a perturbation in the particle orbits, the fraction 77 can be writte n 
in terms of the energy content in the resonant MHD waves (e.g., iBlandford Eichler 1 119871 ). 
7] = (5B I B)~ 2 , of amplitude 5B. We can then use Equation (11) to constrain 6 and the level 
of turbulence. For most cases, we found that k ~ 0.02 which requires that 86° < 9 < 90° and 
6 < r] < 16. This implies a relatively high turbulence level, 

0.25 < -J^ < 0.4. 



4.4. Forward shock and Precursor 

The results indicate that most of the radiation is originated from behind the forward shock 
(see Figure [2|). However, the model could not explain the observed linear profile of the filaments 
(see Figure [5] for an example). The model predicts a sharp decline after the peak, which is not 
observed. 



We speculated that some of the emission may come from a precursor (jEllison et al .11 19941 ). We 
estimated the contribution from the precursor by requiring that the distribution function should be 
continuous across the shock. The precursor would consist of particles that have diffused across the 
shock but remain energetic. Specifically, it evolves in the presence of a magnetic field consistent 
with that of the un-shocked medium surrounding the SNR, here assumed to be 4 times lower than 
the compressed field estimated downstream; however, we assumed that k remains the same. 

Addition of this component substantially improves the predicted profile shape, as shown in 
Figure [5) On the other hand, we found that the inclusion of a precursor hardly affects the spectral 
parameters. More details will be discussed in a future publication (D. Lomiashvili et al. 2010, in 
preparation). 
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5. Discussion 

From the power-law fits to the spectra of the filaments in Cas A, it is seen that the emission 
from the inner regions is consistently softer, by about 10%, than that from the outer regions. This 
seems to be consistent with the effect expected from radiative cooling, since the outer regions have 
had less time since they interacted with the shock. When only synchrotron losses and advection are 
taken into account, however, the predicted difference between the inner and outer photon indices 
is the same in all filaments. The data show that this difference can change from one filament to 
another. 

Also, from a consideration of synchrotron losses, one might expect that the widths of the 
filaments get narrower at higher energies. In fact, if synchrotron cooling and advection were the 
only processes controlling the plasma distribution, the width of a non-thermal filament can be 
estimated as w = V a ^ Y T s , or the distance the particles are advected before radiating away their 
energy. This can also be written as w = V^y/ipB 2 ^, where tp = ar/^mc, with ot the Thompson 
cross section for electrons, B the magnetic field, and 7 the Lorentz factor of the accelerated particle, 
which when assuming emission peaked at the Larmor frequency vl can be written as (v/vl) 1 ^ ■ 
Therefore, an important dependence of the widths on the frequency of the radiation, of the form 
w a zv -1 / 2 , would result. However, the data suggest that no dependence exists. 

These observations seem to point at the existence of additional mechanisms affecting the 
evolution of the plasma distribution and are found to be consistent with the model used. In this 
model, the difference between the photon indices of the inner and outer regions is regulated by 
diffusion, and it is determined mostly by the ratio of advection to diffusion lengths, (, whereas 
it is found that varying A^jf (= l^/R) produces changes mainly in the calculated width of the 
non-thermal filaments without considerably affecting the model spectra. This was also seen when 
carrying out the fits, since the values of x 2 did not change appreciably for a wide range of values 
of A^f, and therefore additional constraints on this parameter were necessary. As a constrain, we 
used the values of A^f that were calculated to match the filamentary widths at half intensity to 
the actual data as initial values for the fits. This parameter ranges from 0.014 to 0.034, while C 
varies from 3.3 to 8.9. 

Due to the role that C plays in the model, it should be possible to correlate the spatial differences 
between the photon indices with the value of the proportionality constant in the diffusion coefficient, 
k. It is seen that Filaments 1, 6, and 7 show the highest values for k (although considerable 
uncertainties were obtained) and that the relative change in photon spectral index from inner to 
outer regions is the lowest for these filaments (although Filament 5 shows a similar change). The 
amount of diffusion can change 1 order of magnitude for the different filaments depending on the 
degree of spatial spectral variation observed. It can be argued that the diffusion of particles tends 
to homogenize the plasma distribution and lower the difference between the inner and outer photon 
indices. This was also seen in the simulations where higher values of k were used. 

The particle spectral indices are found close to 3, although steeper values are also seen (up to 
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4 for Filament 6). This index corresponds to the power-law distribution of electrons resulting from 
shock acceleration processes. After the particles evolve in the magnetic field, one might expect 
to see steeper spectral indices, especially at higher X-ray energies. Such steep distributions might 
have been seen before. For instance, when comparing the 10-32 keV RXTE Proportional Counter 
Array spectrum of Cas A with the predicted synchrotron emission spectrum dominating the band 
from 0.3 keV to 7 keV, the observed excess can be accounted for by an additional contribution of 
non-thermal Bremsstrahlung radiation from a s teep power-law (in dex ~4.1) population of electrons 
interacting with other electrons and with ions (jAllen et al.ll2008l ). 



Besides accounting for these observations, we estimated the magnetic field in each filament 
as well as the level of turbulence. The derived magnetic field is on average below the values 
previously inferred for the NE rim (Filaments 1-3) by assuming that the widt h of the filament 
(take n as 1".5 — 4") is determined by synchrotron losses and advection only (jVink fc Laming 
20031 ). -Bsync ~ 80 — 160^G. The difference might be explained by the fact that our estimates take 



into account these two processes but additionally consider diffusion. The turbulent magnetic field is 
constrained to be 0.25 < ^ < 0.4. Moreover, the ordered field is found to be nearly perpendicular 
to the shock front (86° < 6 < 90°), consistent with an expansion inside a toroidal external field 
produced originally by the progenitor star. The inferred magnetic field in the filaments is still 
higher than that expected from magnetic field amplification of the interstellar field (~ 3/uG), as 
was also pointed out by Vink & Laming. Perhaps the interstellar field surrounding Cas A is higher, 
or the field has been amplified by the high-energy particles near the shock front through nonlinear 
wave growth. 



In some cases, the mag n etic fi eld m ight be much higher in oth e r filam ents of Cas A (jAtovan et al 



2000). iPatnaude &: Fesen I (|2007l ) and lUchivama &: Aharonian I (|2008l ) have observed X-ray vari- 
ability of some of the non-thermal filaments seen in projection near the reverse shock on a time 
scale of a few years. Assuming that this variability is related to synchrotron cooling and DSA, and 
that the diffusion is close to the Bohm limit with the field parallel to the shock front, Uchiyama & 
Aharonian estimated that the field required would be ~ 1 mG. However, there is still debate as to 
whether the observed variability is associated with filaments in t he forward shock, and there are 



other plausible scenarios for explaining it (e.g.. lBykov et al.ll2008l ). 



Models describing the non-equilibrium ion ization behind an SNR forward shock recently devel- 
oped (lEllison et al. 1 120071 : IPatnaude et al. I2OO9I ) show that efficient DSA could increase the ionization 
fraction of some elements. In their model, Patnaude et al. (2009) have found that ionization occurs 
more rapidly and closer to the shock when the particles experience efficient shock acceleration, 
meaning that spatial variations in the intensity of the thermal emission near the shock might be 
used to probe the properties of cosmic rays. We looked at the spectral features and estimated the 
equivalent widths of Si K XIII and S K XV in order to quantify any possible spatial variations in 
their intensity. In Table IU we summarize the results. Unfortunately, the error bars associated with 
the equivalent widths that we observed are too large for us to draw definitive conclusions. 
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Finally, our model assumes that synchrotron radiation is monochromatic at v c (see Equation 
(1)), which is clearly an oversimplification. When we incorporated the full synchrotron spectrum 
into the model, we saw changes in the model parameters. For instance, a fit to the spectrum of 
Filament 5 with the revised model leads to B ~ 70//G and k ~ 0.06, which are different to the 
values shown in Table However, the changes do not qualitatively modify our conclusions. The 
obliquity angle is still close to 90° and the turbulent field is ^ 0.2. The details of the full model 
will be presented in a future publication. 



6. Summary 

We summarize our main results as follows. 



1. Spectral evolution is seen across non-thermal filaments in Cas A, with the spectra of the outer 
regions being harder by about 10% on average. 

2. The widths of the filaments show no significant dependence on photon energy. 

3. To account for the observational results (1 and 2), we needed to include the effects of diffusion. 
If we restrict to Bohm-type diffusion, we could quantify the level of turbulence (0.25 < ^ < 
0.4) as well as the diffusion itself (k ~ 0.02 — 0.15). Moreover, we found that the magnetic 
field is of the order of tens of /iG, varying from filament to filament, and that the field is 
nearly perpendicular to the shock front. 



These results are in overall agreement with models of cosmic ray acceleration in the shocks of 
SNRs. They imply that there is a high level of magnetic turbulence in the non-thermal filaments 
associated with the forward shock of Cas A, as well as magnetic field amplification. Both of these 
conditions are necessary to efficiently accelerate cosmic rays. 

Regarding the shock orientation, our analysis of a sampling of non-thermal filaments, which 
have good azimuthal coverage of the remnant, implies that the obliquities are close to 90°, which 
i s consistent with the ex pansion of Cas A in the wind environment produced by the progenitor 



(IChevalier fc Oishi 112003^ . 



Efficient cosmic ray acceleration in the shock of Cas A would have implications regarding 
the acceleration of protons (and heavier ions), which may interact with cold ambient protons and 
produce neutral pions that would decay into gamma rays, leaving a signature in the spectrum of 
the remnant in the GeV to TeV energy range. This signature could in principle be detected. 



We thank M. Laming, M. Pohl, and S. Reynolds for useful discussions. This research has made 
use of data obtained from the Chandra Data Archive and the Chandra Source Catalog, and software 
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Table 1: Widths of Linear Profiles of Non-Thermal Filaments. 



Filament 0.3 - 2.0 keV 3.0 - 6.0 keV 

>+0.07 
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5 


0.95 
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1.28 
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9 
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+0.13 
-0.11 
+0.09 
-0.07 
+0.26 
-0.20 
+0.03 
-0.02 
+0.07 
-0.06 
+0.08 
-0.06 
+0.17 
-0.06 
+0.24 
-0.18 



1.23 
1.07 



+0.07 
-0.07 
+0.11 
0.10 

1 n3+°- 07 

i - u,3 -0.07 

1 64+ - 26 

-L-^-O.lS 

0.781835 

o.98i8:8 8 7 
i.08i8:8 8 7 
0.9218:8^ 
i-solgff 



6.0 - 1 0.0 ke V 

1+0.41 



1.24] 
1.25 



0.24 
+0.37 



-0.27 

1 65+ 39 
i - DO -0.63 

77+°- 14 

U. I /_ .12 

1 25+ 6 - 1 

1-2618:82 
o.92i8:l 



Note. — Error intervals at the 90% confidence level and widths are in arcseconds. 



Table 2: Best-Fit Spectral Power Laws. 
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Note. — Error intervals at the 90% confidence level 
"Power- law index 

6 Power-law normalization, in units of 10~ 5 photons keV -1 cm" 2 s 
c Reduced chi-squared of the fit 
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Table 3: Key Parameters of the Diffusion- Advection Model. 
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Note. — Error intervals at the 90% confidence level 



Table 4: Equivalent Widths of Si K XIII and S K XV Emission Lines 
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Error intervals at the 90% confidence level 




Fig. 1. — Exposure-corrected image of Cas A in the range from 0.3 keV to 10 keV. The selected non- 
thermal filaments are indicated, along with source and background regions for spectral extraction. 
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Fig. 2. — Linear intensity profile of Filament 5 in the energy range 0.3 keV-10 keV. The top panel 
shows a Chandra image of the filament with the inner and outer regions labeled and the position 
of the shock indicated. 
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Fig. 3. — Linear profiles of the filaments near the peak, for three energy bands (from top to bottom): 
0.3 — 2 keV, 3 — 6 keV, and 6 — 10 keV. The solid lines show the best-fit Gaussian functions. 
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Fig. 4. — Spectra of the inner and outer regions of non-thermal filaments. The solid lines show the 
best power-law fits to the spectra. The residuals of the fits are also shown. 
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Fig. 5. — Fit to the linear profile of Filament 5 (also see Figure 2) obtained by the model described 
in Section 4 (dashed line) and calculated with the addition of a precursor (solid line). 



